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Using 14M '4){2S) events taken with the BES-II detector, Xcj 2{K'^K~) decays are studied. For the four- 
kaon final state, the branching fractions are Z3(xcO,i,2 ^ 2{K+K-)) = (3.48 ± 0.23 ± 0.47) x 10"^, (0.70 ± 0.13 ± 
0.10) X 10~^, and (2.17 ± 0.20 ± 0.31) x 10"^. For the (pK+K' final state, the branching fractions, which are 
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1. Introduction 

Exclusive quarkonium decays provide an im- 
portant laboratory for investigating perturbative 
quantum chromodynamics. Compared with J /ip 
and ip{2S) decays, there is much less knowledge 
on Xcj decays which have parity and charge con- 
jugation PC = ++. Relatively few exclusive de- 
cays of the XcJ have been measured. For the 
Xcj Vector Vector (VV) mode, measurements 
of XcJ (t><l> [1], K*{892fK*{S92f [2], and 
LULo [3] have been recently reported. The search 
for new decay modes and measurements with 
higher precision will help in better understand- 
ing various Xcj decay mechanisms [4,5] and the 
nature of ^Pj cc bound states. 

Furthermore, the decays of XcJ, especially 
Xco £md Xc2, provide a direct window on glue- 
ball dynamics in the 0++ and 2++ channels 
since the hadronic decays may proceed via 
cc ^ gg ^ qqqq. Recently, a paper by 
Zhao [6] points out that the decay branching 
fractions for scalar glueball candidates (/o(1370), 
/o(1500), and /o(1710)) in Xco decays may be 
predicted by a factorization scheme, in which 
some parameters can be fitted with B{XcJ — * 
ujuj,K*{892fK*{m2f,cj)(j)) [6]. The measure- 
ment precision of B{xcj VV) will affect the 
uncertainties of the fitted parameters. Also, these 
fitted parameters will help clarify the role played 
by OZI-rule violation and SU(3) flavor breaking 
in the decays. 

In this analysis, XcJ — * K~^K~K~^K~ is stud- 
ied using tp{2S) radiative decays. The branching 
fractions of XcJ K^K^K^K^ and Xco.2 4><f> 
are measured with higher statistics, and those of 
XcJ decaying to (j)K^K~ are measured for the 
first time. 

2. The BES detector 

The Beijing Spectrometer (BES) is a conven- 
tional solenoidal magnet detector that is de- 
scribed in detail in Ref. [7] ; BESII is the upgraded 
version of the BES detector [8]. A 12-laycr vertex 
chamber (VC) surrounding the beam pipe pro- 
vides trigger and position information. A forty- 
layer main drift chamber (MDC), located radi- 



ally outside the VC, provides trajectory and en- 
ergy loss (dE/dx) information for charged tracks 
over 85% of the total solid angle. The mo- 
menta resolution is (Tp/p = 0.017-^/1 +p'^ (p in 
GeV/c), and the dE/dx resolution for hadron 
tracks is ~ 8%. An array of 48 scintillation coun- 
ters surrounding the MDC measures the time- 
of-flight (TOF) of charged tracks with a reso- 
lution of ~ 200 ps for hadrons. Outside of 
the TOF counters is a 12-radiation- length bar- 
rel shower counter (BSC) composed of gas tubes 
interleaved with lead sheets. This measures the 
energies of electrons and photons over ^ 80% of 
the total solid angle with an energy resolution of 
as/E = 22%/^/e {E in GeV). Outside of the 
solenoidal coil, which provides a 0.4 Tesla mag- 
netic field over the tracking volume, is an iron flux 
return that is instrumented with three double lay- 
ers of counters that identify muons of momenta 
greater than 0.5 GeV/c. 

A GEANT3 based Monte Carlo (MC) program 
with detailed consideration of the detector perfor- 
mance (such as dead electronic channels) is used 
to simulate the BESII detector. The consistency 
between data and Monte Carlo has been care- 
fully checked in many high purity physics chan- 
nels, and the agreement is quite reasonable [9]. 

3. Event selection 

The data sample used for this analysis consists 
of (14.00±0.56) X 10^ ■ip{2S) events [10] collected 
with the BESII detector at the center-of-mass en- 
ergy = M^^^s)- The XcJ ^ K+R-R+R- 
channels are investigated using ip{2S) radiative 
decays to XcJ- Events with four charged tracks 
and one to three photons are selected. Each 
charged track is required to be well fitted by a 
helix and to have a polar angle, 9, within the 
fiducial region | cos 6*1 < 0.8. To ensure tracks 
originate from the interaction region, we require 

Vxy = \Jv[+y^ < 2 cm and \V,\ < 20 cm, 
where V^, Vy, and Vz are the x,y, and z coor- 
dinates of the point of closest approach of each 
charged track to the beam axis. All charged 
tracks must be identified as kaons using the com- 
bined dE/dx and TOF information. 

A neutral cluster is considered to be a photon 
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Figure 1. distribution of 4-C fit for Xco Figure 2. M^|^_ versus M]^|^_ for candidate 

K~^K~K~^K~ . Histogram denotes MC and error events, 
bars denote data. 



candidate if it is located within the BSC fidu- 
cial region (| cos0| < 0.75), the energy deposited 
in the BSC is greater than 50 MeV, the first hit 
appears in the first 6 radiation lengths, the angle 
between the cluster and the nearest charged track 
is more than 15°, and the angle between the di- 
rection of cluster development and the direction 
of the photon emission is less than 40°. 

A four constraint (4-C) kinematic fit under 
the ipi'^S) — > 'yK~^K~K^K~ hypothesis is per- 
formed, and the oi the fit is required to be 
less than 25. For events with two or three pho- 
ton candidates, the combination having the mini- 
mum is chosen. In addition, X^k+k- k+k- ^ 

are required to suppress background contamina- 
tion from i(}{2S) — > 77r+7r~7r+7r~ and '0(25') 
■jTT^TT^ . Figure 1 shows the distribu- 
tion of data and MC for the process of XcO ~^ 
K+K-K+K-. 

With four selected kaons, there are four ways 
to combine oppositely charged kaons, and two 
combinations of M^|^_ m]^|^„ pairs can be 
formed. The combination that has one of its 
Mx+K~ closest to the mass is selected for fur- 
ther analysis. 

Figure 2 shows the distribution of AI^+k- 



versus Mj^+j^- for selected events. There are 
two clear bands near 1.02 GeV/c^ which corre- 
spond to the (j>K^K~ final state. The insert 
in the upper right corner, is the enlarged view 
of the lower left corner, and a clear (j)(f) signal 
can be seen. Figure 3 shows the correspond- 
ing m]^|^_ versus Af]^|^_ distribution for MC 
■0(25) 7Xci, Xci (l)K+K^ and 00, using 
Bixci ^ (I^K+K-) : Bixci ^ 00)=2:1. 

To investigate intermediate resonances in 
[K^K^) final states, the invariant mass distri- 
bution of all four possible K~^K~ combinations 
are plotted in Figure 4(a). Except for the 0, no 
other obvious resonance is seen. 

To test if the above selection criteria will 
cause "fake" signals, "0(25) 7Xci,Xci 

MC events are generated accord- 
ing to phase space. Figure 4(b) shows the 
Mk+k- distribution of these events using the 
same selection as for data. No peak is seen around 
the signal region. 

4. MC simulation 

For each of the channels studied 100,000 MC 
events are generated. The proper angular distri- 
butions for the photons emitted in ^0(25) ^XcJ 
are used [11]. Phase space is used for the Xc.i 
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Figure 3. M'^'_^j^_ versus 
i,{2S) ^ 7Xci, Xci -> 
after event selection. 
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Mk+k-k+k- distribution from this channel and 
that from data, shown in Figure 6, the contribu- 
tion from ip{2S) 'jK^K^K^K^ background 
should be small. No Xcj peaks are seen in Fig- 
ure 5(a) and (b). Thus, the fitted number of Xcj 
signal events is insensitive to the shape of the 
function used to describe the total background. 

We have also considered possible background 
from ^b{2S) ^ 7XcO,XcO ^ /o(980)/o(980) 

. Using the branching ratio for 
Xco ^ /o(980)/o(980) ^ tt+tt-K+K' [13], 
the ratio B(/o(980) ^ K+ K-)/{B{f a{980) 
K+K-) + B{fo{980) TT+TT-)) [13], and MC ef- 
ficiencies for these processes, the estimated num- 
ber of /o(980)/o(980) -> K+K-K+K~ events is 
about 6, which are spread over a relatively wide 
region compared with the width of the (j) [14]- 
Thus, this background in the </> region is negligi- 
bly small. 



2{K^K ) decays (including intermediate states, 
e-g- XcJ 

5. Background study 

Possible backgrounds come from ili{2S) — > 
7Xcj, XcJ Tr+TT^vr+vr^, h+tt~K+K^ , KsKtt, 
and TT+TT-pp] ii){2S) -> n°K+K-K+K-; and 
phase space 'ip{2S) -> -fK+K-R+R-. 100,000 
MC events are generated for each of the first 
four background channels where the angular dis- 
tribution for the radiative photon is generated 
the same as for the signal channels while the 
XcJ decays are generated according to phase 
space. The contamination from these back- 
grounds to each signal channel is less than 1.0%, 
which is negligible. 100,000 MC events are 
also generated for each of the last two back- 
ground channels according to phase space. For 
i){2S) Tr° R+R'R+R', Figure 5(a) shows the 
Mk+k- K+K- distribution in the region 3.2-3.7 
GeV/c^. Using the branching fraction measured 
by CLEOc [12], the number of events from this 
channel is expected to be 23 ± 7. For ip{2S) 
-fR+R-R+R-, the 

K-K+K- distribution 
is shown in Figure 5(b). This branching fraction 
is currently unavailable. However, comparing the 



6. Mass spectrum fit 

6.1. XcJ ^ 2{K+K-) 

Figure 6(a) shows the high mass (> 3.2 
GeV/c^) 2{R~^R~) invariant mass distribu- 
tion using all candidate events in — > 
lXcJ,XcJ 2{R+R^). Clear XcJ signals can 
be seen in this Figure. A fit with Breit-Wigner 
functions convoluted with Gaussian resolution 
functions (about 15 MeV/c^ for Xcj signals in 
all channels studied) yields the number of Xcj 
events: N^^^ = 278 ± 18, N^^^ = 54 ± 10, and 
Ny^^^ = 160 ± 14. MC simulation gives detection 
efficiencies of 5.9%, 6.3%, and 5.8% for J = 0, 1, 
and 2, respectively. 

6.2. XoJ (t)K+K- 

The XcJ (j)R^R^ events are clustered 
as two bands around the mass in Figure 2. 
The (t)R+R- region is defined by m]^^+^_ e 

(1.08,2.4) GeV/c2 and M]^f^_ e (1.00,1.04) 
GeV/c^ (i, j = 1, 2). Figure 6(b) shows the 
^^ct>K+K- distribution for these events, and clear 
XcJ signals can be seen. MC simulation gives de- 
tection efficiencies of 5.9%, 6.2%, and 5.6% for 
J=0, 1, and 2, respectively. A fit with Breit- 
Wigner functions convoluted with Gaussian reso- 
lution functions yields the number of Xcj events: 



iV^.o = 53.5 ± 8.3, iV^^i = 19.2 ± 5.6, and 
iV^„2 = 56.3 ± 8.2. 

In order to determine the contribution from 
non-resonant 'ip{2S) ■yXcJ,XcJ 2{K+K-) 
events in the (j)K~^K~ region, we analyze the non- 
resonant 2{K^K^) region, defined by m'^\^_ S 
(1.1,2.4) GeV/c^ (i = 1,2) in Figure 2. The num- 
ber of events in this region is TV"'' = 238. Fig- 
ure 7(a) shows the M^+k-k+k- distribution for 
the non-resonant 2{K'^ K^) region events, where 
three clean XcJ peaks arc seen with very little 
background. A fit with Breit-Wigner functions 
convoluted with Gaussian resolution functions 
yields: N^.. = 144.1 ± 12.7, N^^. = 21.5 ± 5.5, 
and iV^- =39.4 ±6.9. 

Figure 7(b) shows the </> signal for the (j)K^K~ 
region events. The fit yields 143 ± 14 signal 
events and Nbg = 25 background events from 
non-resonant 2{K~^K~) events within a 3a win- 
dow (1.00 - 1.04 GeV/c^) in the (p signal re- 
gion. We calculate the non-resonant 2{K~^K~) 
contribution to N^^g, -^xci; ^^'^ -^Xc2 assum- 
ing that the proportion of XcOj Xci, and Xc2 
events is the same for non-resonant 2{K^K~) 
events in the non-resonant and the (j)K~^K~ re- 
gions. For Xco, the non-resonant contribution is 
N^^g = N^^lg ■ Nbg/N^"- = 15.1 ± 1.3. Similarly, 
we obtain TV" , = 2.3 ± 0.6, and AT" = 4.1 ± 0.7. 
Therefore, the number of (f)K'^K~ events, after 
subtracting the non-resonant 2{K'^K~) events, 
is N^^„ = 38.4 ± 8.4, N^^^ = 16.9 ± 5.6, and 
iVxc2 = 52.2 ± 8.2. 

6.3. XoJ (fxp 

The signal region for XcJ — > 4'4> events is a 40 
McV/c^x 40 McV/c^ square aromid the <j) mass 
in Figure 2. Figure 8 shows the M^^ distribution 
for these 00 events, and clear Xco,2 signals can be 
seen. MC simulation gives detciction efficiencies 
of 9.0% and 8.6% for J=0 and 2, respectively. A 
fit yields the number of events for Xcj- TVxco = 
27.8 ± 5.8 and iV^cs = 42.7 ± 7.1. 

The two bands near the mass in Figure 2, 
used to extract the Xcj — * (pK~^K~ signal in Sec- 
tion 6.2, are taken as the sideband region for the 
00 events. They include both the (j)K^K~ events 
and non-resonant K'^K~K'^K~ events. From 
MC simulation the event distributions in the two 



bands are nearly uniform. The number of nor- 
malized sideband events in the 00 signal region 
are 

respectively, where the factor f=33 is the ratio of 
the sideband area to the 00 signal region area. 
Thus we obtain the number of events in XcJ 
N^^g = 26.2 ± 5.8 and iVxc2 = 41.0 ± 7.1. 

7. Systematic error 

The systematic error in these branching frac- 
tion measurements includes the uncertainties 
caused by wire resolution, particle ID, photon ef- 
ficiency, and the number of '0(25') events. 

The systematic error caused by MDC track- 
ing and the kinematic fit are estimated by using 
simulations with different MDC wire resolutions 
[9]. For particle ID, the combined information of 
dE/dx and TOF is used. An error of 2% is as- 
signed for each charged track [16] and each pho- 
ton [9] . The errors introduced by branching frac- 
tions of intermediate states are taken from the 
Particle Data Group (PDG) [15]. 

The total systematic errors, determined by the 
sum of all sources added in quadrature, are listed 
in Table 1. The uncertainty from ,B(0 K^K~) 
contributes once in the systematic error estima- 
tion for XcJ — * 4'K^K~ and twice in 00, while 
it docs not contribute in Xc.j 2{K^K~). For 
the uncertainties caused by wire resolution, there 
are some slight differences for the different decay 
channels. 

8. Results 

For XcJ 2{K^K~) (including intermediate 
states), the branching fractions are calculated us- 
ing 

N 

B{xcJ - 2{K+K-)) = — -fl^ — , 

N^(2S) ■ B{ip{2S) jxcj) ■ e 



Table 1 

Systematic error (%). In the wire resolution row, the numbers from left to right correspond to ^^{28) 
2{K+K-), (j)K+K-, and (jxj). 



Source 


XcO 


Xcl 


Xc2 


Wire resolution 


8.9, 9.8, 10.0 


9.3, 9.9 


9.7, 9.6, 10.1 


Particle ID 


8 


8 


8 


Photon efficiency 


2 


2 


2 


Background shape 


negligible 


negligible 


negligible 


Number of ij;{2S) 


4 


4 


4 


B(^(25) ^ 7Xcj) 


4.3 


4.6 


4.9 




K+K-) 


1.2 


1.2 


1.2 


Total 


XcJ - 2{K+K-) 

XcJ ^ <PK+K- 
XcJ (p<P 


13.5 

14.1 
14.3 


13.9 
14.3 


14.3 

14.2 
14.5 



where, the average detection efficiency e is given 
by 

e = e2(K+K-) + 

^4>K+K- + -Tl- ^<t»t>- 

Similarly, we can calculate the branching frac- 
tions for XcJ — * (f)K~^K~, (l)(p with corresponding 
efficiency expressions. Tabic 2 lists our measure- 
ment results, together with the PDG values. 

In summary, the decays of XcJ — * 2{K~^K~) 
arc studied, and the corresponding branching 
fractions including intermediate states are given. 
The decay XcJ — * (pK~^K~ is observed for the 
first time. The branching fractions for Xcj 
2{K~^K~) and XcJ 4"t' are measured with 
higher precision; Table 2 lists the comparison of 
the measured branching fractions between BESII 
and the PDG. Our measurement for XcJ ~^ 4"Py 
together with the two measurements of XcJ — * 
and i?*(892)"/f* (892)0, ^ie helpful in under- 
standing the nature of XcJ states. 
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Figure 4. Mx+k- distribution of all four possible 
K^K~ combinations, (b) Mj^+k- distribution 
for MC V(25) ^ 7Xci,Xci ^ K+K-R+R-. 
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Figure 5. Background shape obtained from MC 
simulation, (a) ij}{2S) -k^R+R-R+R'. (b) 
V'(25) ^ -fR+R-R+R-. 




Figure 6. Breit-Wigner fit to Xcj signals (a) with 
all candidate events and (b) with the (j)K^K~ 
region events. 
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Figure 8. Fit to XcJ signals in (jxp final state. 



Figure 7. (a) Breit-Wigner fit to XcJ using non- 
resonant 2{K^K^) region events, (b) Fit to ip 
signal with (j)K^ K~ region events. 



